Research in Astronomy and Astrophysics, 24:085014 (12pp), 2024 August 
© 2024. National Astronomical Observatories, CAS and IOP Publishing Ltd. Printed in China and the U.K. 


https: / /doi.org/ 10.1088 /1674-4527 /ad621d 


CrossMark 


Pulsations of Three Rapidly Oscillating Ap Stars TIC 96315731, TIC 


72392575, and TIC 318007796 


Hai-Jian Zhong’, Dong-Xiang Shen!, Chun-Hua Zhu', He-Lei Liu!, Su-Fen Guo!, and Guo-Liang Li? 
School of Physical Science and Technology, Xinjiang University, Urumqi 830046, China; chunhuazhu @sina.cn 


Received 2024 April 27; revised 2024 June 6; accepted 2024 June 25; published 2024 August 14 


Abstract 


We analyze the frequencies of three known roAp stars, TIC 96315731, TIC 72392575, and TIC 318007796, using 
the light curves from the Transiting Exoplanet Survey Satellite. For TIC 96315731, the rotational and pulsational 
frequencies are 0.1498360 day ' and 165.2609 day” ', respectively. In the case of TIC 72392575, the rotational 
frequency is 0.25551day '. We detect a quintuplet of pulsation frequencies with a center frequency of 
135.9233 day ', along with two signals within the second pair of rotational sidelobes of the quintuplet separated 
by the rotation frequency. These two signals may correspond to the frequencies of a dipole mode. In TIC 
318007796, the rotational and pulsational frequencies are 0.2475021 day ', 192.73995day ', and 
196.33065 dag `. respectively. Based on the oblique pulsator model, we calculate the rotation inclination (i) 
and magnetic obliquity angle (8) for the stars, which provide the geometry of the pulsation modes. Combining the 
phases of the frequency quintuplets, the pulsation amplitude and phase modulation curves, and the results of 
spherical harmonic decomposition, we conclude that the pulsation modes of frequency quintuplets in TIC 
96315731, TIC 72392575, and TIC 318007796 correspond to distorted dipole mode, distorted quadrupole mode, 
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and distorted dipole mode, respectively. 
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1. Introduction 


Chemically peculiar stars with spectral type A, known as Ap 
stars, typically exhibit magnetic field with intensities ranging 
from 300 G to 34 kG (Babcock 1960; Auriére et al. 2007; Elkin 
et al. 2010). The intense magnetic field suppresses convection 
in the stellar interior, which leads to an interaction of 
gravitational settling and radiative levitation in the atmosphere 
of Ap stars (Michaud 1970). The intricate interaction results in 
their surfaces manifesting an overabundance of elements such 
as Pr, Nd, Eu, Gd, and Tb, which can reach values up to a 
million times that of the solar value (Ryabchikova et al. 2004; 
Michaud et al. 2015). The overabundances of rare earth 
elements form chemical spots near their magnetic poles 
(Ryabchikova et al. 2007), which remain stable over periods 
ranging from days to decades (Renson & Manfroid 2009). 
From the light curves modulated by the spots, we can 
determine the rotation period of the stars. 

Rapidly oscillating Ap (roAp) stars are a rare subset within 
the Ap stars (Kurtz 1982), located at the lower boundary of the 
classical instability strip and intersecting the main sequence. 
These stars exhibit high-overtone (n > 15), low-degree (J < 3) 
non-radial pressure mode (p-mode) pulsations with periods 
ranging from 4.7 to 25.8 minutes (e.g., Cunha et al. 2019; 
Holdsworth et al. 2021). In the Johnson B band, the semi- 


amplitudes of these pulsations can reach up to 18 mmag 
(Holdsworth 2015). In the majority of roAp stars, pulsations are 
believed to be driven by the «-mechanism within the HI 
ionization zone (e.g., Dziembowski & Goode 1996; Gautschy 
et al. 1998; Balmforth et al. 2001; Cunha 2002). However, for 
some roAp stars with pulsation frequencies exceeding the 
theoretical acoustic cut-off frequency limit, their pulsations are 
thought to be driven by turbulent pressure in the envelope 
region where convection is not suppressed (Cunha et al. 2013). 

Kurtz (1982) introduced the oblique pulsator model for roAp 
stars, which was further developed by incorporating the 
influence of both rotation and magnetic fields (e.g., Dziem- 
bowski & Goode 1985; Shibahashi & Takata 1993; Takata & 
Shibahashi 1994, 1995; Bigot & Dziembowski 2002; Saio & 
Gautschy 2004; Bigot & Kurtz 2011). In this model, the 
pulsation axis is misaligned with the rotation axis, typically 
aligning with the magnetic axis. As an roAp star rotates (with a 
non-zero rotation inclination angle), the observed pulsation 
amplitude and phase undergo modulation. This modulation 
causes the pulsation frequency to split into 2/+ 1 frequency 
components, each separated by an interval equivalent to the 
rotation frequency (Kurtz 1982). This not only provides 
geometric information about the pulsation modes but also 
facilitates pattern recognition. This peculiarity of roAp stars 
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makes them excellent subjects for studying the interaction 
between magnetic fields and stellar pulsations. 

After Kurtz (1978) discovered the first roAp star, the search 
for roAp stars has been ongoing (e.g., Martinez et al. 1991; 
Balona et al. 2011; Holdsworth et al. 2014, 2024; Hey et al. 
2019; Holdsworth 2021). In these search efforts, ground-based 
telescopes face some disadvantageous aspects in detecting low- 
brightness and low-amplitude roAp stars due to Earth’s rotation 
and atmospheric effects. However, space telescopes such as 
Kepler (Koch et al. 2010) and the Transiting Exoplanet Survey 
Satellite (TESS; Ricker et al. 2015), with continuous and high- 
precision observation capabilities, can effectively overcome 
these disadvantages. Especially TESS, with its all-sky coverage 
and high-precision short-cadence photometric observation 
capabilities, has propelled the exploration of roAp stars. The 
latest observations from TESS indicate that the confirmed 
count of roAp stars has now reached 118 (Holdsworth et al. 
2024). Considering the extremely limited number of known 
roAp stars, a reevaluation of the pulsations exhibited by these 
acknowledged roAp stars using TESS’s latest observational 
data will aid in revealing potential pulsation signals that may 
exist but have yet to be detected, thus enabling a more 
comprehensive understanding of the pulsation phenomena 
associated with roAp stars. 

In this work, we report the results of frequency analysis, 
pulsation mode identification, and the geometry of pulsation 
modes for three roAp stars: TIC 96315731, TIC 72392575, and 
TIC 318007796. In Section 2, we introduce the three identified 
roAp stars, along with the TESS observations. A detailed 
analysis of frequency characteristics is provided in Section 3, 
while Section 4 presents pulsation modes and their geometry. 
The characteristics of pulsation amplitude and phase modula- 
tion curves are discussed in Section 5, and Section 6 presents 
the spherical harmonic decomposition results of the amplitude 
spectra. Finally, we present the conclusions for the stars in 
Section 7. 


2. The roAp Stars and TESS Observations 
2.1. Three Known roAp Stars 


TIC 96315731 (HD 51203) is classified as an oC CND star with 
a spectral type of A2p SrEuCr (Houk 1982; Renson et al. 1991). 
According to the report by Buder et al. (2018), the effective 
temperature of the star is 7098 + 105 K. Holdsworth et al. (2021) 
identified distinct rotational and pulsational features in the light 
curves from TESS sectors 6 and 7, officially designating the star 
as anew roAp star. They determined the stellar rotation period as 
6.6713 + 0.0007 days, consistent with the value reported by 
Bernhard et al. (2015) of 6.6752 + 0.0008 days. Additionally, 
Holdsworth et al. (2021) determined the frequency of the central 
peak in the frequency quintuplet is 165.258 + 0.001 day ' 
(1.91271 + 0.000 01 mHz). Chojnowski et al. (2019) reported a 
mean magnetic field modulus of 7.9 + 0.5 kG for TIC 96315731 
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by analyzing magnetically split absorption lines in the H-band 
spectra obtained from the Apache Point Observatory Galactic 
Evolution Experiment (APOGEE) survey. 

TIC 72392575 (HD 225578) exhibits an ambiguity in its 
spectral classification. Initially, it was classified as ApSr: by 
Cowley & Cowley (1964), but later Renson et al. (1991) 
determined it as A5 Sr. Paunzen et al. (2012) reported the 
effective temperature of the star to be 8290+ 60K. Paunzen 
et al. (2012) used photometric data from the 1m Austrian- 
Croatian Telescope with a Bessell B filter system over 
130.9 minutes and did not detect any pulsation signals, with 
the corresponding maximum amplitude in the Fourier spectrum 
found to be 1.6 mmag. Based on the TESS light curves, Balona 
(2022) initially identified this star as an roAp star and derived its 
rotation period as 3.922 days. Recently, Holdsworth et al. (2024) 
detected a pulsation signal in the TESS sector 14 light curve at a 
frequency of 135.944 + 0.003 days (1.57343 + 0.000 03 mHz), 
with no evidence of pulsation frequency rotational sidelobes. 
They also determined the rotation period of the star as 
3.9016 + 0.0008 days. 

TIC 318007796 (HD 190290) is classified as an o CND star 
with a spectral type of A2 EuSrCr (Kazarovets et al. 1993; 
Polosukhina et al. 1999). The Strémgren—Crawford indices 
measured by Joshi et al. (2016) for the star are as follows: 
b—y=0.289, m,=0.293, c,;=0.466 and H= 2.796, 
giving 6m,;=—0.091 and Ac =—0.306. According to the 
report by Shen et al. (2023), the effective temperature of 
the star is 7079 + 163 K. Martinez & Kurtz (1990) obtained 
high-speed photometric data for 8.85 hr at the South African 
Astronomical Observatory by the 1.0m Elizabeth telescope 
with a Johnson B filter. They identified rapid oscillations at 
frequencies 196.13 day~! (2.27 mHz) and 192.67 day~! (2.23 mHz), 
revealing a frequency separation of 40 «Hz between the two 
signals. Furthermore, after pre-whitening these frequencies, 
the amplitude spectrum hinted at additional unresolved 
frequencies. van Heerden et al. (2012) conducted differential 
photometric observations and derived a rotation period of 
4.03 + 0.13 days but did not mention pulsation frequencies. 
Recently, Holdsworth et al. (2024) detected pulsation 
signals with frequencies of 192.73997 + 0.00013 days 
(2.2307867 + 0.000 001 5 mHz) and 196.33370 + 0.00002 days 
(2.2723808 + 0.000 000 2 mHz) in TIC 318007796, using 
light curves from TESS sectors 27 and 39. The rotation period 
of the star was determined to be 4.09070 + 0.00002 days. 
They emphasized that the mode with the lower frequency 
corresponds to a distorted dipole mode, while the mode at the 
other frequency represents a distorted quadrupole mode. This 
suggests a large frequency separation of approximately 
83 uHz. Furthermore, Hubrig et al. (2004) and Bagnulo 
et al. (2006, 2015) conducted measurements on the mean 
longitudinal magnetic field of TIC 318007796. Based on these 
measurements, Rustem et al. (2023) reported an averaged 
quadratic effective magnetic field of (Be) = 3104 + 182 G. 
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Table 1 
Detailed Information on TESS Light Curves for TIC 96315731, TIC 72392575, and TIC 318007796 


Name Sector BJD Time Span (days) Number of Points Exposure Time (s) 
TIC 96315731 6 2458468.27520-2458490.04775 21.77 14,639 120 
TIC 96315731 7 245849 1.63525-2458516.08916 24.45 16,084 120 
TIC 96315731 33 2459201.73606-2459227.57521 25.84 17,458 120 
TIC 96315731 34 2459229.19326-2459254.06781 24.87 16,689 120 
TIC 72392575 54 2459769.90234—2459796.04423 26.14 17,417 120 
TIC 72392575 55 2459797.10394—2459824.26598 27.16 18,882 120 
TIC 318007796 27 2459036.28194—2459060.64411 24.36 16,781 120 
TIC 318007796 39 2459361.77426-2459389.72065 27.95 19,337 120 


Note. The first column gives the stellar names, the second column lists the light curve-associated sectors, and the third to fifth columns provide the start and end times 
of the light curves, the duration, and the count of points, respectively, with the sixth column specifying the exposure time. 


2.2. TESS Observations and Data Reduction 


Our targets are observed in different sectors by TESS. In 
detail, TIC 96315731 is observed in sectors 6, 7, 33, and 34; 
TIC 72392575 in sectors 54 and 55; and TIC 318007796 in 
sectors 27 and 39. The detailed information on the light curves 
for these sectors is shown in Table 1. We employ the Python 
package Lightkurve (Lightkurve Collaboration et al. 2018) to 
download 2 minute cadence light curves of the stars processed 
by the Science Processing Operations Center from the Mikulski 
Archive for Space Telescopes server. 

The PDC_SAP fluxes of the light curves are used in our 
analysis. The Lightkurve package is also employed to remove 
outliers greater than 5g from the light curves, as well as other 
points with quality issues, such as sudden changes in 
background measurements and impulsive outliers. The light 
curves for each sector are then divided by their median flux to 
achieve flux normalization, and the fluxes are subsequently 
converted to relative magnitudes. For TIC 96315731, we 
combine the light curves of sectors 6, 7, 33, and 34, resulting in 
a total of 64,780 data points over an observational period of 
96.94 days. Similarly, for TIC 72392575, sectors 54 and 55 are 
combined, yielding 36,299 data points observed over 
53.30 days. For TIC 318007796, sectors 27 and 39 are 
combined, comprising 36,118 data points observed over 
52.31 days. Figure 1 displays the combined light curves of 
the stars. It is evident that the combined light curves of the stars 
display distinct rotational features, which are controlled by the 
spots on their surfaces. 


3. Frequency Analysis 


We perform a Discrete Fourier Transform (DFT; Kurtz 
1985) on the combined light curves of TIC 96315731, TIC 
72392575, and TIC 318007796 to obtain their rotation and 
pulsation frequency information. The Nyquist frequency 
ua = 1/(2At), where At is the sampling time interval of 
TESS observations, is the highest frequency without under- 
sampling at the given sampling rate. For TESS observations, 


the Nyquist frequency is approximately 359 day '. Conse- 
quently, the frequency range of the amplitude spectrum spans 0 
—359day '. When the frequency separation between two 
frequencies exceeds the frequency resolution, defined as 
Jes = 1.5/T, they become distinguishable (Loumos E 
Deeming 1978). Here, T is the time span of the combined 
light curve. For TIC 96315731, TIC 72392575, and TIC 
318007796, their respective frequency resolutions are 
0.01547 day—', 0.0281 day~', and 0.02868 dan '. 

Figure 2 illustrates the amplitude spectra of TIC 96315731, 
TIC 72392575, and TIC 318007796. Panels (a) display the 
amplitude spectra covering the frequency range of 
0-359 day |. Panels (b) exhibit signals related to the rotational 
frequency, while panels (c) show signals associated with the 
pulsation frequency. The residual amplitude spectra after pre- 
whitening the pulsation frequency are presented in panels (d). 
The red dashed lines represent 4.6 times the background signal, 
while the gray dashed lines signify the extracted frequencies. 


3.1. Rotation Frequency Analysis 


During the analysis of the rotation frequency, we apply the 
DFT to extract both the rotation frequency and its harmonics 
within the frequency range of 0—2 day’. Following this, linear 
least squares fitting is employed to determine the amplitude and 
phase values associated with the identified frequencies. Lastly, 
optimization of the frequency, amplitude, and phase results is 
conducted using the nonlinear least squares fitting technique 
(same as the technique in Kurtz 1982). The results of the 
nonlinear least squares fitting for TIC 96315731, TIC 
72392575, and TIC 318007796 are presented in Table 2. 


3.1.1. TIC 96315731 


For TIC 96315731, the amplitude spectrum in the low-frequency 
region is shown in Figure 2(b). The results of the rotation frequency 
and its first harmonic are presented in Table 2. The derived rotation 
frequency of the star is Vot = 0.1498360 + 0.0000005 den ` 
(Bot = 6.67396 + 0.000 02 days). It is noteworthy that the 
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Figure 1. Top panel: Light curves of TIC 96315731 with four different sectors; Middle panel: Light curves of TIC 72392575 with two different sectors; Bottom panel: 


Light curves of TIC 318007796 with two different sectors. 


formal error reported here is calculated using the method 
proposed by Montgomery & O’Donoghue (1999), which 
represents the lower limit of the error. Our rotation period is 
consistent with the findings reported by Bernhard et al. 
(2015), who measured a rotation period of 6.6752 + 
0.0008 days, and by Holdsworth et al. (2021), who recorded 
a rotation period of 6.6713 + 0.0007 days. 


3.1.2. TIC 72392575 


For TIC 72392575, the amplitude spectrum in the low- 
frequency region is depicted in Figure 2(b). Table 2 presents the 
results of the rotation frequency and its first harmonic. The 
derived rotation frequency of the star is Vot = 0.25551 + 
0.00002 day"! (Ro = 3.9138 + 0.0004 days). Our result is 
consistent with the rotation period of 3.9016 + 0.0008 days 
determined by Holdsworth et al. (2024). Furthermore, our 


analysis benefits from a longer overall observational period, 
which results in reducing uncertainty in the derived rotational 
period. 


3.1.3. TIC 318007796 


For TIC 318007796, Figure 2(b) displays the amplitude 
spectrum in the low-frequency region. The results of the 
rotation frequency and its first harmonic are shown in 
Table 2. The derived rotation frequency of vot = 0.2475021 + 
0.0000004 day~' (Pot = 4.040370 + 0.000 007 days). Within 
the uncertainties, our rotation period is consistent with the one 
derived by van Heerden et al. (2012), which is 4.03 + 0.13 days. 
However, it is smaller than the rotation period of 
4.09070 + 0.00002 days obtained by Holdsworth et al. (2024). 
This discrepancy may mainly come from the different research 
methods. 
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Figure 2. The amplitude spectra of TIC 96315731, TIC 72392575, and TIC 318007796. (a) The amplitude spectra within the frequency range from 0 to 359 day™'. (b) 
Amplitude spectra within the frequency range of 0—2 day '. (c) Amplitude spectra of the region containing pulsation frequency after pre-whitening of low-frequency 
signals. (d) The residual amplitude spectra after pre-whitening the pulsation frequency. The red dashed lines indicate 4.6 times the white noise, while the gray dashed 
lines signify the extracted frequencies. 
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Table 2 
Nonlinear Least Squares Fitting Results for All Extracted Frequencies in TIC 96315731, TIC 72392575, and TIC 318007796 
Star Label Frequency Amplitude Phase 
(day~') (mmag) (rad) 
TIC 96315731 iot 0.1498360 + 0.0000005 10.40 + 0.01 3.154 + 0.001 
Jet 0.2996231 + 0.0000014 3.43 + 0.01 2.190 + 0.003 
V — Jet 164.95976 + 0.00005 0.063 + 0.007 —1.26 + 0.12 
V — Vrot 165.10972 + 0.00002 0.161 + 0.007 1.98 + 0.05 
DN 165.26090 + 0.00003 0.108 + 0.008 —1.93 + 0.07 
V + Yor 165.41075 + 0.00002 0.139 + 0.007 —1.15 + 0.05 
V + Vrot 165.56060 + 0.00005 0.061 + 0.007 2.69 + 0.13 
TIC 72392575 Vrot 0.25551 + 0.00002 4.19 + 0.01 2.289 + 0.002 
rot 0.51166 + 0.00002 3.91 + 0.01 2.277 + 0.003 
n 135.4916 + 0.0023 0.022 + 0.005 —1.41 + 0.23 
Vz 135.7454 + 0.0009 0.056 + 0.005 —0.90 + 0.09 
V3 — dier 135.4089 + 0.0019 0.026 + 0.005 2.36 + 0.19 
V3 — Veo 135.6659 + 0.0010 0.053 + 0.005 2.41 +0.10 
V3 135.9233 + 0.0003 0.190 + 0.005 2.60 + 0.03 
V3 + Vrot 136.1785 + 0.0010 0.052 + 0.005 2.41 +0.10 
v3 + 2brot 136.4358 + 0.0024 0.021 0.005 2.41 + 0.24 
TIC 318007796 Veot 0.2475021 + 0.0000004 23.208 + 0.009 2.7989 + 0.0004 
Jet 0.4950086 + 0.0000020 4.282 + 0.009 1.2453 + 0.0020 
EIN 0.9900158 + 0.0000070 1.219 + 0.009 2.9984 + 0.0072 
Vi — Vrot 192.48973 + 0.00006 0.115 0.007 1.21 + 0.06 
v 192.73995 + 0.00013 0.050 + 0.007 —0.94 + 0.14 
Vi + Vrot 192.98758 + 0.00007 0.094 + 0.007 —1.65 + 0.07 
Vy = Wrot 195.83561 + 0.00011 0.061 + 0.007 2.66 £0.11 
V2 — Vrot 196.08629 + 0.00008 0.085 + 0.007 1.33 + 0.08 
Vz 196.33065 + 0.00002 0.310 + 0.007 —3.07 + 0.02 
V2 + Vrot 196.57818 + 0.00006 0.106 + 0.007 —1.81 + 0.06 
Vy + D 196.82260 + 0.00006 0.111 + 0.007 0.33 + 0.06 


Note. The time zero-points for fitting to TIC 96315731, TIC 72392575, and TIC 318007796 are BJD 2458861.12908, BJD 2459797.32875, and BJD 2459213.62661, 
respectively. The table columns list star names, frequency labels, frequencies, amplitudes, and phases, in that order. 


3.2. Pulsation Frequency Analysis 


Before determining the pulsation frequencies, we perform 
high-pass filtering on the combined light curves of the stars to 
remove the effects of low-frequency signals such as stellar 
rotation and instrumental long-term trends. This ensures the 
accuracy of analyzing high-frequency pulsation signals. High- 
pass filtering is achieved by continuously pre-whitening peaks 
within the frequency range of 0-20 day! until their amplitude 
levels are comparable to the noise level in the high-frequency 
region. Prior to high-pass filtering, careful examination is 
conducted to ascertain the absence of astrophysical signals 
within the high-frequency region. We utilize the batch mode of 
the Period04 software (Lenz & Breger 2005) to execute the 
high-pass filtering. Similar to the extraction of rotation 
frequencies, we perform DFT on the high-pass filtered 
combined light curves to obtain pulsation frequencies for each 
star. Amplitude and phase values of the pulsation frequencies 
are then obtained through linear least squares fitting, followed 
by the application of nonlinear least squares fitting to optimize 


the frequencies, amplitudes, and phases (same with the 
technique as in Kurtz 1982). 


3.2.1. TIC 96315731 


The amplitude spectrum after high-pass filtering, as depicted 
in Figure 2(c), reveals that TIC 96315731 exhibits a quintuplet 


of pulsation frequencies centered at v= 165.2609 day 
Additionally, rotational sidelobes are observed at v 4 


=1 


E Vrot and 


LZ, The results for pulsation frequencies are presented in 


Table 2. After pre-whitening the quintuplet of frequencies, a 
slight power excess is observed in the residual amplitude 
spectrum, as shown in Figure 2(d). The power excess suggests 
the presence of unresolved signals in the amplitude spectrum. 


3.2.2. TIC 72392575 


From Figure 2(c), the amplitude spectrum of TIC 72392575 
after high-pass filtering indicates the presence of three pulsation 
frequencies: vı = 135.4916 day `. Vy = 135.7454 day ', 
and a quintuplet of pulsation frequencies centered at 
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v = 135.9233 day |. The peak of the second pair of rotational 
sidelobes associated with 13 is lower than 4.6 times the 
background noise level (but remains above 4 times the 
background noise level). Given that their frequency spacing 
from the first pair of rotational sidelobes closely approximates 
the rotational frequency, we have incorporated them into our 
frequency detection results. The results for pulsation 
frequencies are presented in Table 2. It is worth noting that 
the frequency separation v — vı = 0.2538 day _' is very close 
to the rotation frequency Vrot = 0.25551 dan `. This suggests 
that vı and 12 may be part of the frequency triplet of a dipole 
mode. After pre-whitening the pulsation frequencies, a slight 
power excess is evident in the residual amplitude spectrum, as 
shown in Figure 2(d). The power excess indicates the presence 
of unresolved signals remaining in the amplitude spectrum. 


3.2.3. TIC 318007796 


After high-pass filtering, the amplitude spectrum reveals the 
existence of two pulsation frequencies for TIC 318007796: a 
triplet centered around a frequency vı = 192.73995 day! and a 
quintuplet centered at a higher frequency vz = 196.33065 day `. 
as depicted in Figure 2(c). The detailed results of the pulsation 
frequencies are presented in Table 2. Within the uncertainty 
range, our results are consistent with those reported by 
Holdsworth et al. (2024). Upon pre-whitening the pulsation 
frequencies, slight power excesses are observed in the residual 
amplitude spectrum within the frequency ranges of 192.0- 
192.5 day~' and 195.0-196.0 day ', as shown in Figure 2(d). 
These power excesses indicate the presence of unresolved signals 
at these locations. 

It is noteworthy that our investigation of TIC 96315731, TIC 
72392575, and TIC 318007796 has exclusively unveiled p- 
mode pulsations indicative of the stellar outer layers, with no 
detection of g-modes sensitive to the core region (Aerts et al. 
2010). Consequently, the task of ascertaining the complete 
rotational profile, spanning from the surface to the core, of 
these particular roAp stars remains a challenge within the scope 
of current asteroseismic analyses. 


4. Pulsation Mode and Geometry 


For roAp stars, the frequency triplet is commonly associated 
with the dipole mode, while the frequency quintuplet is 
typically linked to the quadrupole mode. However, some roAp 
stars, such as HD 6532 (Kurtz et al. 1996) and HR 3831 (Kurtz 
et al. 1990), respectively exhibit the dipole mode as the 
pulsation modes of the frequency quintuplet and frequency 
septuplet. This is due to the influence of the magnetic field, 
which can increase the degeneracy of the pulsation mode, 
allowing dipole mode to exhibit five or more fine structures 
(Saio & Gautschy 2004). Therefore, solely relying on the 
number of rotational sidelobes in the pulsation frequency does 
not allow us to determine the pulsation mode of roAp stars. In 
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this paper, we assume that the pulsation mode associated with 
the quintuplet of pulsation frequencies is a quadrupole mode. 
According to the oblique pulsator model, the rotational 
sidelobes of pure dipole or quadrupole mode are precisely 
split by the rotation frequency, with equal phase for all 
components of the multiplet at the maximum pulsation 
amplitude. To test this, we employ the rotation frequency 
fixed rotational sidelobes corresponding to the pulsation 
frequency. It is ensured that the interval between these 
rotational sidelobes precisely equates to the rotation frequency. 
Subsequently, a suitable time zero-point is chosen such that the 
phases of the first pair of rotational sidelobes for the pulsation 
frequency are equal, Le, (V — Mo) = OU + Mot). Finally, 
linear least squares fitting is employed to the combined light 
curves to obtain phase values for all components. The results of 
the linear least squares fitting for TIC 96315731, TIC 
72392575, and TIC 318007796 are presented in Table 3. 

According to the oblique pulsator model, there are the 
following relationships between the amplitude ratio of the 
quintuplet of pulsation frequency and the stellar rotation 
inclination (i) and magnetic obliquity angle (8) for the 
quadrupole mode (Kurtz et al. 1990): 


AR +AQ ` 3 sin? d sin? i De 
AP (3cos? 6 — 1)(3cos?i — 1)’ 
A? + A2 12 sin 8 sin i cos 8 cos i (2) 


A  GBcos?B — I Gcos?i — 1)” 


From Equations (1) and (2), we can derive the standard 
constraints for the oblique pulsator 


2 4 4Q 
Ary + ASS 


tani tan 3 = 4———___ 
(2) 2)? 
AÌ 4 AM 


(3) 


where A®, AP), and AB represent the amplitudes of the 
central peak and the first and second pair of rotational sidelobes 
for the frequency quintuplet, respectively. The tani tan 8 
expression provides constraints on the stellar rotation 
inclination (i) and magnetic obliquity angle (8), but these 
constraints are limited. We cannot directly determine the values 
of i and Ø from these constraints. Simultaneously, solving 
Equations (1) and (2) allows us to obtain the values of i and £. 
These provide the geometry of the pulsation mode. 


4.1. TIC 96315731 


For TIC 96315731, the phases of the frequency quintuplet 
are unequal, as shown in Table 3. This indicates that the mode 
corresponding to the frequency quintuplet is distorted. 
Additionally, the amplitudes of each pair of rotational sidelobes 
are unequal. This is attributed to the effect of the Coriolis force 
(Bigot & Dziembowski 2002). By substituting the amplitudes 
from Table 3 into Equations (1), (2), and (3), we determine that 
tani tan 8 = 1.63 + 0.14. Here, the stellar rotation inclination 
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Table 3 
Linear Least Squares Fitting Results for Pulsation Frequencies of TIC 96315731, TIC 72392575, and TIC 318007796, with Frequency Spacing Equal to the Rotation 
Frequency 
Star Label Frequency Amplitude Phase 
(day~!) (mmag) (rad) 
TIC 96315731 V — rot 164.96122 0.061 + 0.007 1.86 + 0.12 
V — Vrot 165.11106 0.160 + 0.007 —1.16 + 0.05 
v 165.26090 0.109 + 0.007 —1.93 + 0.07 
V + Vrot 165.41074 0.140 + 0.007 —1.16 + 0.05 
V + 2r%0¢ 165.56058 0.061 + 0.007 2.71 + 0.12 
TIC 72392575 Vy 135.4915 0.022 + 0.005 —1.43 + 0.23 
V2 135.7448 0.054 + 0.005 —0.91 + 0.09 
V3 — Wot 135.4123 0.024 + 0.005 2.37 + 0.21 
V3 — Vrot 135.6678 0.052 + 0.005 2.42 + 0.10 
V3 135.9233 0.191 + 0.005 2.60 + 0.03 
V3 + Vrot 136.1788 0.052 + 0.005 2.42 + 0.10 
V3 + Wot 136.4343 0.021 + 0.005 2.38 + 0.24 
TIC 318007796 V1 — Vrot 192.49245 0.108 + 0.007 —1.91 + 0.06 
Vy 192.73995 0.049 + 0.007 —0.95 + 0.14 
Vi + Vrot 192.98745 0.093 + 0.007 —1.67 + 0.07 
Vz — Wot 195.83565 0.062 + 0.007 2.65 + 0.11 
V2 — Vrot 196.08315 0.083 + 0.007 —1.81 + 0.08 
V2 196.33065 0.310 + 0.007 —3.07 + 0.02 
Vz + Vrot 196.57815 0.107 + 0.007 —1.81 + 0.06 
Vz + Vrot 196.82565 0.111 + 0.007 —2.78 + 0.06 


Note. The time zero-points for fitting to TIC 96315731, TIC 72392575, and TIC 318007796 are BJD 2458861.12908, BJD 2459797.32875, and BJD 2459213.62661, 


respectively. The table is organized with columns for star names, frequency labels, frequencies, amplitudes, and phases. 


angle is i= 68°52 + 0°92 and the magnetic obliquity angle is 
B= 32°62 + 1°29, or vice versa. 


4.2. TIC 72392575 


For TIC 72392575, the phases of the frequency quintuplet 
are not equal but very close, as depicted in Table 3. This 
indicates a slight distortion in the mode corresponding to the 
frequency quintuplet. Similar to TIC 96315731, based on the 
amplitudes of the quintuplet from Table 3 and Equations (1), 
(2), and (3), we determine that tani tan 8 = 1.73 + 0.30. Here, 
the stellar rotation inclination angle is i= 78°08 + 1°15 and 
the magnetic obliquity angle is G=20°07+ 1°50, or 
vice versa. 


4.3. TIC 318007796 


For TIC 318007796, the phases of the frequency quintuplet 
are unequal, as shown in Table 3. This implies that the mode 
corresponding to the frequency quintuplet is distorted. Similar 
to TIC 96315731, according to the amplitudes of the 
quintuplet from Table 3 and Equations (1), (2), and (3), we 
derive that tani tan 8 = 3.64 + 0.28. Here, the stellar rotation 
inclination angle is i=81°43+0°51, and the magnetic 
obliquity angle is G= 28°77 + 0°66, or vice versa. Within 
the uncertainty range, our result of ont on is consistent 


with the finding of Holdsworth et al. (2024), where they 
reported tani tan 6 = 3.51 + 0.27. 


5. Pulsation Amplitude and Phase Modulation 


According to the oblique pulsator model, the pulsation nodes 
for a pure quadrupole mode are located at a colatitude of 
+54°7, while for a pure dipole mode, the pulsation nodes 
coincide with the equator. If i+ 8 > 54°7 (or 90°, for a pure 
dipole mode) the pulsation nodes will cross the line of sight, 
resulting in alternating visibility of two pulsation poles. This 
causes the pulsation amplitude to be zero, accompanied by a 
phase reversal of 7-rad. By analyzing the pulsation amplitude 
and phase modulation curves, we can determine whether there 
is distortion in the pulsation mode. In addition, the variations in 
the rotational light are governed by the abundance spots on the 
stellar surface, while the changes in pulsation amplitude are 
controlled by the angle between the pulsation axis and the line 
of sight. Therefore, the phase difference between the maximum 
values of the rotational light curve and the pulsation amplitude 
indicates the relative positioning of the abundance spots and 
the pulsation poles. 

To test this, we divide the combined light curves into 
segments of specified length based on the pulsation period. 
Then, we apply linear least squares fitting to each segment to 
obtain the amplitude and phase values corresponding to the 
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Figure 3. Phase folded light curves (top panels), the modulation curves of the pulsation amplitude (middle panels), and phase (bottom panels) as a function of rotation 
phase for TIC 96315731, TIC 72392575, and TIC 318007796. The bin size of the phase folded light curves is 0.005. Larger phase errors occur at smaller pulsation 
amplitudes, hence points with phase errors greater than 1o have been excluded. The red curves represent the fitted curves of pulsation amplitude and phase modulation 
using the method proposed by Kurtz (1992). The time zero-points for TIC 96315731, TIC 72392575, and TIC 318007796 are BJD 2458861.12908, BJD 


2459797.32875, and BJD 2459213.62661, respectively. 


pulsation frequency. These amplitude and phase data are 
employed to construct pulsation and phase modulation curves 
as a function of the rotation phase. 


5.1. TIC 96315731 


For TIC 96315731, we divide the combined light curve into 
129 segments, with each segment having a length of 120 
pulsation cycles, corresponding to a duration of 17.43 hr. Then, 
we apply linear least squares fitting to each segment, obtaining 
the amplitude and phase values corresponding to the pulsation 
frequency of 165.2609 day’. Figure 3 shows the rotational 
light curve, as well as the pulsation amplitude and phase 
modulation curves as a function of the rotation phase. At 
rotation phase 0, the phase of maximum pulsation amplitude 
closely approximates that of the maximum rotational light. This 
implies that the abundance spots are close to the pulsation poles 
(assuming the abundance spots are near the magnetic poles). At 
this time, the angle between the pulsation axis and the line of 
sight is i— G= 35°90. Following a half-cycle rotation, the 
secondary maximum of the pulsation amplitude occurs, with 
the angle between the pulsation axis and the line of sight being 
i+ @=101°14. Furthermore, near the rotational phase 0.25, 
the pulsation amplitude reaches its minimum, accompanied by 
a phase reversal of ad in the pulsation phase. These 
characteristics are consistent with the expectations of the 
oblique pulsar model. For the quadrupole mode (or dipole 


mode), the minimum pulsation amplitude is non-zero, 
suggesting a distortion in the pulsation mode. 


5.2. TIC 72392575 


For TIC 72392575, we divide the combined light curve into 
129 segments, each with a length of 100 pulsation cycles, 
corresponding to approximately 17.77 hr. Then, we perform a 
linear least squares fitting on each segment, obtaining the 
amplitude and phase values corresponding to the pulsation 
frequency of 135.9233 day’. Figure 3 displays the rotational 
light curve, as well as the modulation curve of the pulsation 
amplitude and phase as a function of rotation phase for the star. 
At rotation phase 0, the phase of the maximum pulsation 
amplitude is not consistent with the maximum of the rotation 
light, which indicates that the abundance spot is not aligned 
with the pulsation poles. This phenomenon is quite common in 
roAp stars, as evidenced by the studies on the roAp star HD 
86181 by Shi et al. (2021) and on HD 12098 by Kurtz et al. 
(2024). At this moment, the angle between the pulsation axis 
and the line of sight is i— 8 = 58°01, with the upper limb of 
the pulsation node appearing on the line of sight. At rotation 
phase 0.25, the pulsation amplitude reaches its minimum 
(which is not zero), with a slight variation in the pulsation 
phase but not a phase reversal of z-rad. This implies that the 
quadrupole mode is distorted. At rotation phase 0.5, the 
pulsation amplitude reaches its second maximum, with the 
angle between the pulsation axis and the line of sight being 
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Table 4 
The Spherical Harmonic Decomposition Results of Frequency Quintuples for TIC 96315731, TIC 72392575, and TIC 318007796 
Star l AQ A® AQ? AM Af} o 
(mmag) (mmag) (mmag) (mmag) (mmag) (rad) 
TIC 96315731 2 0.061 0.150 —0.109 0.150 0.061 1.860 
1 e 0.309 0.368 0.289 e —1.219 
0 Wi 0.400 e e 2.084 
TIC 72392575 2 0.024 0.054 —0.191 0.05 0.021 2.37 
1 e 0.003 0.004 0.003 e —1.737 
0 os Wi 0.381 Wi e 2.477 
TIC 318007796 2 0.062 0.082 —0.310 0.108 0.111 2.650 
1 e 0.130 0.082 0.170 e —1.158 
0 Wi 0.647 Wi e 2.828 


Note. The time zero-points for TIC 96315731, TIC 72392575, and TIC 318007796 are BJD 2458861.12908, BJD 2459797.32875, and BJD 2459213.62661, 
respectively. The rotation inclination (i) and magnetic obliquity (8) angles are i = 68°52 and 8 = 32°62 for TIC 96315731, i = 78°08 and 8 = 20°07 for TIC 
72392575, and i = 81°43 and 8 = 28°77 for TIC 318007796. The table is organized into columns listing star names, spherical harmonic degrees (/) for different 
pulsation modes, amplitudes for azimuthal orders ranging from m = —2 to m = +2, and phases for each spherical harmonic degree. 


i+ 8=98°15. The pulsation phase between 0.25 and 0.75 in 6. Spherical Harmonic Decomposition 
the pulsation phase modulation curve appears to be more 
scattered compared to other positions. This phenomenon can be 
attributed to the inverse relationship between the uncertainty of 
pulsation phase and pulsation amplitude (Montgomery & 
O’Donoghue 1999). We note that the pulsation amplitude and 
phase modulation curves of the star closely resemble those of 
two other roAp stars, KIC 10685175 (Shi et al. 2020) and HD 
86181 (Shi et al. 2021), pulsating in distorted quadrupole 
modes. This implies that our assumption may be correct. 


To further determine the contribution of the amplitude of 
pulsation modes associated with different spherical harmonic 
components (J=0, 1, 2, 3) to the amplitude spectra of these 
stars, we employ the axisymmetric spherical harmonic 
decomposition method for roAp star pulsation mode proposed 
by Kurtz (1992). This method utilizes a linear combination of 
axisymmetric (m = 0)/ = 0, 1, 2, and 3 spherical harmonics to 
describe the amplitude spectrum of the pulsation mode. For the 
frequency quintuplet, it is sufficient to characterize the 
amplitude spectrum by employing a linear combination of 
5.3. TIC 318007796 spherical harmonic components with <2. The spherical 
harmonic decomposition results for the frequency quintuplets 
in TIC 96315731, TIC 72392575, and TIC 318007796 are 
presented in Table 4. 


Similar to the previous two stars, we divide the combined 
light curve of TIC 318007796 into 198 segments, each with a 
length of 50 pulsation periods, corresponding to a duration of 
6.11 hr for each segment. We then apply linear least squares 


fitting to ages segment at the pulsation frequency of 6.1. TIC 96315731 

196.33065 day ~ to obtain the corresponding amplitude and 

phase values. The rotational light curve, along with the Our results for TIC 96315731 indicate that the dipole (/ = 1) 
pulsation amplitude and phase modulation curves are shown mode contributes 0.966 mmag to the pulsation amplitude. 
in Figure 3. At rotation phase 0, the approximate equal phase Additionally, the radial (7 = 0) mode contributes 0.400 mmag, 
between the maximum pulsation amplitude and the maximum while the quadrupole (/ = 2) mode contributes 0.313 mmag. 
of the rotating light indicates that the abundance spots are close Among all the modes, the dipole mode dominates the 
to the pulsation poles. In this case, the angle between the contribution of the amplitude, indicating that the pulsation 
pulsation axis and the line of sight is i— G= 52°66. Around mode of the star is a distorted dipole mode rather than a 
rotation phase 0.25, the pulsation amplitude reaches its distorted quadrupole mode. It is worth noting that our result 
minimum (which is not zero). Additionally, there is a variation differs from that reported by Holdsworth et al. (2021) for TIC 
in the pulsation phase, but it is less than ad This phase 96315731. In their work, they leaned toward a distorted 
change indicates that the pulsation mode is distorted. After quadrupole mode based on the analysis of pulsation amplitude 
rotation phase 0.25, the pulsation amplitude reaches its second and phase modulation curves. According to the spherical 
maximum, which is approximately equal to the maximum harmonic decomposition results presented in Table 4, we 
pulsation amplitude. The angle between the pulsation axis and model the pulsation amplitude and phase modulation curve for 
the line of sight becomes i + 8 = 110°20. TIC 96315731, illustrated by the red curves in Figure 3. 
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6.2. TIC 72392575 


For TIC 72392575, our results reveal that the contributions 
of the radial (1 =0) mode, dipole (J = 1) mode, and 
quadrupole (J = 2) mode to the pulsation amplitude are 
0.381 mmag, 0.01 mmag, and —0.042 mmag, respectively, at 
the maximum of pulsation amplitude. The radial mode 
dominates contribution of the amplitude, followed by the 
quadrupole and dipole modes. Because the pulsation modes of 
the roAp stars are non-radial, the contribution of the dipole 
mode to the amplitude can be neglected compared to the 
quadrupole mode. Therefore, our hypothesis is valid, indicating 
that the pulsation mode of the star is a distorted quadrupole 
mode. We model the pulsation amplitude and phase modulation 
curves for the star using the spherical harmonic decomposition 
results presented in Table 4. The results are shown by the red 
curves in Figure 3. 


6.3. TIC 318007796 


According to the spherical harmonic decomposition results, 
at the maximum pulsation amplitude, the radial (J = 0) mode, 
dipole (J = 1) mode, and quadrupole (/ = 2) mode contribute 
0.647 mmag, 0.382 mmag, and 0.053 mmag to the pulsation 
amplitude, respectively. After excluding the radial mode, the 
contribution of the dipole mode to the pulsation amplitude is 
greater compared to the quadrupole mode. Therefore, we 
conclude that the pulsation mode corresponding to the 
frequency quintuplet is a distorted dipole mode. Utilizing the 
spherical harmonic decomposition results from Table 4, we 
model the pulsation amplitude and phase modulation curves for 
the star, which are shown in Figure 3 as red curves. 


7. Conclusions 


TIC 96315731 is an roAp star with a quintuplet of pulsation 
frequencies. We determine its rotation frequency to be Uu zz 
0.1498360 + 0.0000005 day~' (Po. = 6.67396 + 0.000 02 days) 
and its pulsation frequency to be v= 165.2609 day_' through 
careful analysis of the combined light curve. Then, we determine 
the rotation inclination angle (i = 68°52) and the magnetic 
obliquity angle (8 = 32° 62) through the amplitude ratio of the 
quintuplet of pulsation frequencies, which provides initial 
constraints on the stellar geometry. 

TIC 72392575 is also an roAp star with a quintuplet of 
pulsation frequencies. We determine its rotation frequency to be 
Lu = 0.25551 + 0.00002 day”! (Po = 3.9138 + 0.0004 days). 
The pulsation frequencies of the star display intriguing 
characteristics. In the second pair of rotational sidelobes 
associated with the pulsation frequency v3 = 135.9233 day ', 
two frequencies are observed, with a separation approximately 
equal to the rotation frequency. These signals may correspond to 
the frequencies originating from a dipole mode. According to the 
amplitude ratio of the quintuplet of pulsation frequencies, we 
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determine the rotational inclination angle (i = 78° 08) and the 
magnetic obliquity angle (8 = 20° 07). 

TIC 318007796 is an roAp star with two pulsation modes. 
We derive its rotation frequency as v,.=0.2475021 + 
0.0000004 das ` (Bo: = 4.040370 + 0.000 007 days), and the 
central peaks of the triplet and quintuplet of pulsation 
frequencies correspond to frequencies vı = 192.73995 day ` 
and v= 196.33065 day ', respectively. From the amplitude 
ratio of the quintuplet of the pulsation frequencies 12, we obtain 
the rotation inclination angle (i = 81°43) and the magnetic 
obliquity angle (8 = 28°77). 

We examine the phases of the frequency quintuplets from 
the linear least squares fitting results, along with the pulsation 
amplitude and phase modulation curves of TIC 96315731, TIC 
72392575, and TIC 318007796. The results suggest the 
presence of distortion in the pulsation modes. Integrating these 
findings with the results of spherical harmonic decomposition, 
we confirm that the pulsation modes of the frequency 
quintuplets corresponding to TIC 96315731, TIC 72392575, 
and TIC 318007796 are distorted dipole modes, distorted 
quadrupole modes, and distorted dipole modes, respectively. 
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